A possible black hole binary in Mkn 501 
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Abstract. A simple binary model for Mkn 501 is considered under the assumption 
that the TeV and X-ray periodicity of ~ 23 days, observed during the flaring state 
in 1997, may be basically interpreted as a doppler-shifted flux modulation due to 
the orbital motion of the relativistic jet. For the typical jet properties inferred from 
emission models, we find an intrinsic orbital period of (6 — 14) yrs and a centre-of-mass 
1/-^ ' distance of (2 — 3.5) x 10^^ cm, the binary thus being a very close system. If the binary 

separation corresponds to that at which gravitational radiation becomes dominant, 
one may obtain a maximum allowed primary mass of ~ 10^ Mq and a corresponding 
maximum allowed secondary mass in the range of (1 — 9) x IO^A/q assuming that 
■ gas accretion occurs on around the salpeter time scale. Interestingly such values agree 

, with the black hole masses expected from merger scenarios. 

o 

^ : INTRODUCTION 

I: 

I in the universe. For, in hierarchical clustering scenarios, giant elliptical galaxies 

(e.g. the host galaxy of Mkn 501) are thought to be products of mergers between 
spiral galaxies. Since the brightest galaxies generally seem to contain massive black 
holes (BH) in their nuclei, a BBHS should be a natural consequence, e.g. [2], [6]. 
In the specific case of the gamma-ray blazar Mkn 501 at 2; = 0.034, the existence 
of such a putative BBHS has been repeatedly invoked as a possible explanation 
of the observed misalignment between pc and kpc-jet [3] or as a reason for the 
special evolution of its X-ray SED [10]. The periodicity of Pobs — 23 days, recently 
observed in the X-ray and 7-ray band of Mkn 501 ( [4], [7]; see Fig. 1) may add 
another piece of evidence on such an internal configuration. Here we show, that 
this periodicity may be interpreted as being basically caused by the orbital motion 
of the relativistic jet emerging from the less massive BH in a close binary system 
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FIGURE 1. Fourier analysis of the TeV and X-ray data from April to July 1997 using a nor- 
malized Lomb periodogram. The periodicity of 23 days is clearly identifiable. Figure courtesy of 
D. Kranich. 



MODELLING 

Assuming that the periodicity observed in the flaring state in Mkn 501 is asso- 
ciated with a relativistically moving feature in the jet (e.g. blob, shock, knot) the 
flux modulation by Doppler boosting is given by the simple formula 



^obs 



where a is the spectral flux index. For a source region following a helical path (see 
Fig. 2) the Doppler factor S{t), and thus the observed flux, is a periodical function 
of time and may be written as 



5{t) 



^1 - (t;2 + nlR^)/c^ 



fj. cos 2 



QkR sin i sin t)/c^ 



where Q^. is the keplerian orbital frequency, R the centre-of-mass distance, i the 
jet inclination angle, c the velocity of light and Vz the outflow velocity. Due to 
the relativistic motion of the source region the observed period appears strongly 
shortened 

-Pobs = (1 + -2^) (1 COS i) Pk . 

c 

Thus, for an observed period of Pobs = 23 days, typical inclination angles i ~ 1 /Th 
and bulk Lorentz factors Ft = 10 — 15, e.g. [5,9], we have an intrinsic orbital period 
-Pfe = (6 — 14) yrs . Using the deflnition of the keplerian orbital frequency, one may 
also derive an equation for the required binary mass ratio [8] 
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FIGURE 2. Simple sketch of a helical jet path caused by the orbital motion in a BBHS. 



with m, M the mass of the smaller /larger BH, respectively, G the gravitational 
constant and / ~ 8 the observed TeV flux ratio between maximum and minimum. 
Applying an additional constraint by assuming that the current binary separation 
d roughly corresponds to the separation at which gravitational radiation becomes 
dominant, we may set an upper limit on the allowed binary masses (note that for a 
more precise estimate the gravitational time scale Tgrav in [8] should be multiplied 
by a factor of 4 and accordingly, the rhs of Eq. (10) by 0.25°'^^, thus enhancing 
the allowed upper limit for m by roughly a factor of 3). These upper limits are 
calculated for the relevant inclination angles and two different spectral indices in 
Tab. 1, where Pk denotes the intrinsic orbital period, ti = Tgrav/4 the gravitational 
lifetime and Pp the period of gravitomagnetic jet precession. The calculation has 
been performed assuming a.) that the dynamical timescale for (Eddington) gas 
accretion is of the order of the salpeter time scale, i.e tgas ~ ts — 5 x 10^ yrs (note 
that for tgas > ts the upper limit for m can be reduced) and b.) that the spectral 
index a = 1.2 (in bracket: a = 1.7). The above mentioned constraint may be 
related to the key aspect that BL Lacs are less luminous radio galaxies with rather 
little signs of cosmological evolution [1]. They may thus correspond to the final 
evolutionary stage of BBHSs with decreased mass accretion rates [10]. 



The results indicate that a BBHS in Mkn 501 may not be improbable since 1.) it 
might represent a natural explanation for the observed periodicity; 2.) the putative 
binary black hole interpretation is supported by galaxy merger evolution; 3.) the 
derived masses are in line with such expectations from merger scenarios, i.e. a 




CONCLUSIONS 



TABLE 1. Upper limit for the maximum binary masses. 



i = l/7b 



1/10 



1/15 



m [10^ Mq] 

M [10^ Mq] 
d [10^6 cm] 

Pk [yrs] 



0.13 

1.11 
2.58 
6.10 
1.57 
0.69 



(0.18) 

(0.92) 
(2.48) 



0.81 

1.44 
5.40 
13.7 
2.04 
5.88 



(0.98) 

(1.31) 
(5.43) 



n [10^ yrs] 
Pp[10^yrs] 



(1.30) 
(0.59) 



(1.86) 
(6.05) 



secondary mass of m ~ 10^ Mq appear to be in agreement with the concept that 
the galaxy swallowed in the merger process was a (minor) spiral galaxy; 4.) its 
expected separation is similar to that suggested in a simple unification scenario 
for BL Lacs and FSRQs [10] ; 5.) its precessional period of ~ 10^ yrs would agree 
with the driving frequency found by Conway & Wrobel [3] in order to explain the 
misalignment of the radio jet on pc and kpc-scale. 

We have also checked for consistency with several observations (e.g. VLBI and 
X-ray (RXTE) data) and discussed the imphcations elsewhere [8]. To be sure, the 

presented model is rather simplified because small changes, for example, in the 
maximum electron Lorcntz factor or the magnetic field along the trajectory of the 
emission region may further add to flux variations. However, the general signature 
according to which a rotating jet in a BBHS may cause a periodical flux modulation 
remains robust. Generally, if the binary hypothesis is correct, the observable period 
in Mkn 501 should remain similar during different outburst phases unless there is 
evidence for a change in the general jet properties (i.e. in the bulk Lorentz factor 
or inclination angle). 
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